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CN . Accurate reconstruction of the mass of a resonance decaying to a pair of 

r leptons is challenging because of the presence of multiple neutrinos from 
^ I T decays. The existing methods rely on either a partially reconstructed 

I ■ mass, which has a broad spectrum that reduces sensitivity, or the coUinear 

approximation, which is applicable only to the relatively small fraction of 
^ , events. We describe a new technique, which provides an accurate mass 

reconstruction of the original resonance and does not suffer from the limita- 
^ I tions of the collinear approximation. The major improvement comes from 

\^ ■ replacing assumptions of the collinear approximation by a requirement that 

OO . mutual orientations of the neutrinos and other decay products are consis- 

^ . tent with the mass and decay kinematics of a r lepton. This is achieved 

^ ! by minimizing a likelihood function defined in the kinematically allowed 

phase space region. In this paper we describe the technique and illustrate 
its performance using Z/7*— >rr and if— j-rr events simulated with the re- 
alistic detector resolution. The method is also tested on a clean sample of 
data Z/7*— )-rr events collected by the CDF experiment at the Tevatron. 
We expect that this new technique will allow for a major improvement in 



X 

ci I searches for the Higgs boson at both the LHC and the Tevatron. 



1. Introduction 

Invariant mass reconstruction is commonly used in experimental searches 
for new physics, such as for the Higgs or Z' bosons, as well as in mea- 
surements of properties of known resonances. This technique is relatively 
straightforward for e+e~, /i"*"/!", or di-jet final states. The accuracy of mass 
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reconstruction in these channels is dominated by the detector resolution 
for lepton or jet momenta. The sensitivity of "mass bump-hunting" analy- 
ses depends critically on how narrow the signal invariant mass distribution 
is compared to the (usually broad) distributions in background processes. 
Unfortunately, this simple strategy is much less effective in searches for 
resonances decaying to a pair of r leptons because the r lepton energy asso- 
ciated with neutrinos escapes detection, and only visible products (leptons 
in the case of leptonic r decays or jets in the case of hadronic r decays) are 
observed in the detector. 

Each r lepton decay involves one or two neutrinos, depending on the 
final state: hadronic (r— + hadrons) or leptonic [r^Vr + li^u where l=e 
or /i). In pp or pp collisions, the full energy of neutrinos cannot be deter- 
mined. Instead, one can only reconstruct a transverse energy imbalance in 
the calorimeter (or missing transverse energy, ^), which is representative 
of the total transverse momentum of all neutrinos in the event. Therefore, 
when two or more neutrinos are produced in the same event, their individual 
transverse momenta and directions cannot be reconstructed. The situation 
in decays of heavy resonances into two r leptons is even more complex. In 
these events, the two r's are often produced "back-to-back" and the missing 
momentum associated with their neutrinos partially cancels out. As a re- 
sult, the invariant mass of a resonance cannot be directly reconstructed from 
the ^ and visible decay products of r leptons. Various techniques exist 
to partially reconstruct the mass of resonances in rr final states. However, 
the reconstructed mass distributions for signal processes are rather broad 
(with long tails and typical core resolutions on the order of ~20%), which 
makes it difficult to separate them from the background and considerably 
reduces the signal significance. This poses a major challenge for the Higgs 
boson searches in the H^tt channel, one of the most important channels 
for discovering a low-mass Higgs boson at the LHC [H, 0], whether in the 
context of the Standard Model or beyond (for example, in supersymmetric 
models). Another challenge in searching for a low- mass Higgs boson in the 
TT channel is the large and irreducible background from Z/^*^tt events. 
This is because the Z/'-j* background is several orders of magnitude larger 
than any expected Higgs signal, and its broad partially reconstructed mass 
distribution completely dominates the signal region (for example, see Fig. [1] 
or reference [sl). Therefore, a major improvement in rr invariant mass, 
Mt-ti reconstruction techniques is needed in order to significantly enhance 
the sensitivity of H^tt searches at the Tevatron and LHC experiments. 
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In this paper, we propose a new method, which substantially improves 
the accuracy of the rr invariant mass reconstruction. We expect it will lead 
to a major improvement in the sensitivity of the Higgs boson searches in the 
iJ— T-rr channel at the Tevatron and LHC. In the next section, we briefly 
review currently used methods. Section [3] describes the new technique and 
illustrates its performance using a Monte Carlo simulation with the realistic 
detector resolution. In Sec. HI we report the results of tests on a clean sample 
of data Z/7*— )-rr events collected by the CDF experiment at the Tevatron. 
Finally, we conclude in Sec. [51 




Figure 1: Example of the transverse mass (left plot) defined as an invariant mass of ^ 
and visible r decay products, and the fully reconstructed mass (right plot) using the 
collinear approximation for three event samples: inclusive Z/7*— >rr and gg^H^rr 
with M^f=115 and 130 GeV/c^. Results are obtained for the fully hadronic tt decay 
mode. Events are simulated with the realistic detector resolution (discussed in Sec. I3.3p . 
All distributions are normalized to the unit area. 



2. Review of the Commonly Used Techniques for tt Mass Recon- 
struction 

The two methods frequently used at hadron colliders either rely on re- 
constructing a partial invariant mass or use the collinear approximation. In 
this section, we review these techniques and discuss their advantages and 
shortcomings. 

2.1. The Transverse Mass Method 

Neutrinos from the r decays escape detection and make it impossible 
to determine the 4-momenta of r leptons and thus Mtt- Therefore, one of 
the simplest and frequently used methods relies on a partial (or reduced) 
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invariant mass reconstruction. Examples include either the invariant mass 
of visible decay products of the two r leptons or the invariant mass of the 
visible decay products and in the event, transverse mass. The latter is 
defined as follows: 

M^(rvisi, rvis2, Pt) = m^isi + "^visa + 

where ^=(^^, 0) and Pvisi,2 5 ''^visi,2 are the momenta vectors and in- 
variant masses of the visible r decay products. The transverse mass provides 
a somewhat better separation from the QCD multi-jet backgrounds with 
fake T signatures, and it is often preferred in data analyses. The advantage 
of this technique is that the partial mass can be defined for all signal events, 
thus preserving the statistical power of the available data. However, ignor- 
ing or not fully accounting for the neutrino momenta biases and broadens 
the reconstructed M^-r distributions, and therefore leads to a significantly 
reduced sensitivity in searches and measurements. This problem is particu- 
larly prominent in the low-mass H^tt search, where the signal cannot be 
separated from the much larger and very broad Z^tt background. This 
is illustrated in the left plot in Fig. [H which shows the transverse mass 
^(''"visi, Tvis2 7 distribution. 

2.2. Collinear Approximation Technique 

The collinear approximation is another frequently used technique [ll, [ij. 
This method was first proposed in reference [3] to reconstruct the invariant 
mass in rr decays of a Higgs boson produced in association with an energetic 
jet. It is based on two important assumptions: that the neutrinos from 
each T decay are nearly collinear with the corresponding visible r decay 
products (i.e., 0;^~0vis and O^c^O^is)', and that the ^ in the event is due 
only to neutrinos. In this case, the total invisible momentum carried away 
by neutrinos in each r decay can be estimated by solving two equations: 

= Pmisi sin cos 0visi + Pmis2 6'vis2 COS 0vis2 

fky = Pmisi sin sin (p^-.^^ + p^,^^ sin O^i^^ sin (p^-,^^ , (2) 

where pr^ and pTy are the x- and y-components of the pr vector, Pmisi 
and Pmis2 are the combined invisible momenta (there can be two i/'s in a r 
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decay) of each r decay, and ^visi 2 4>visi 2 ^^e polar and azimuthal 
angles of the visible products of each r decay. Then, the invariant mass 
of the rr-system can be calculated as MT-T-=m^is/ \/^\^2 , where m^is is the 
invariant mass of visible r decay products, and xi,2=Pvisi.2/(Pvisi.2 +Pmisi,2) 
are momentum fractions carried away by visible r decay products. Despite 
offering the great advantage of a fully reconstructed rr mass (Mtt) instead 
of a partial visible mass, the coUinear approximation still has significant 
shortcomings. The technique gives a reasonable mass resolution only for 
the small fraction of events where the rr system is boosted, i.e., produced 
in association with a large jet, and the visible r decay products are not 
back-to-back in the plane transverse to the beam line. The last requirement 
is needed, because the system of Eqs. [2]becomes degenerate if 0visi=0vis2 
and solutions Pniisi,2~sin~"'^(0visi — 0vis2) diverge as |0visi — 0vis2l~^^- Un- 
fortunately, the majority of if— J-rr events are produced with r leptons in 
nearly the back-to-back topology. Therefore, this technique is applicable 
only to a relatively small fraction of rr events. The collinear approxima- 
tion is also very sensitive to the ^ resolution and tends to over-estimate 
the rr mass, leading to long tails in the reconstructed mass distribution 
(see right plot in Fig. [T]). This effect is especially undesirable for low- mass 
Higgs boson searches, where the tails of a much larger Z— j-rr background 
completely overwhelm the expected Higgs peak region. 

3. The Missing Mass Calculator Technique 

The new technique proposed in this paper, the Missing Mass Calculator 
(MMC) method, allows for a complete reconstruction of event kinematics 
in the rr final states with significantly improved invariant mass and neu- 
trino momentum resolutions. The MMC technique does not suffer from 
the limitations of the collinear approximation described in the previous sec- 
tion and can be applied to all rr event topologies without sacrificing the 
reconstructed mass resolution. 

3.1. The Concept and Method Description 

To facilitate the description of the method, we begin with assuming a 
perfect detector resolution and that there are no other neutrinos in rr events 
except for those from the r lepton decays. Under these assumptions, full 
reconstruction of the event topology requires solving for 6 to 8 unknowns: 
X-, y-, and z-components of the invisible momentum carried away by neu- 
trino(s) for each of the two r leptons in the event, and, if one or both r's 
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decay leptonically, the invariant mass of the neutrinos from each leptonic r 
decay. However, there are only 4 equations connecting these unknowns: 

fh:^ = Pmisi sin 6'misi cos 0misi + Pmis2 siu 6'mis2 COS 0mis2 
fky = Pmisi sin sin (f)^;^^ + Pinis2 siu 6'mis2 siu 0mis2 

= mil,, + mli^^ + 2^/p2~Tm2~yp2~^^^^ 

Sj^visj^PiYiisj^ COS A^jiffjj 

^r2 = "^mis2 + "^vis2 + 2y^Pvis2 + "^vis2 \/Pmis2 + "^mis2 

2pvis2Pmis2 COS A^^mj (3) 

where and pTy are the x- and ^/-components of the pr vector, Pvisi,2! 
'^visi 2 5 ^visi 2 5 4'visi 2 are the momenta, invariant masses, polar and azimuthal 
angles of the visible r decay products, and M^=1.777 GeV/c^ is the r lepton 
invariant mass. The rest of the variables constitute the "unknowns" which 
are the combined invisible ( "missing" ) momenta Pmisi 2 carried away by the 
neutrino (or neutrinos) for each of the two decaying r leptons and the 
invariant mass of the neutrino(s) in the r decay, mmig^ j. Finally, S9yrni,2 is 
the angle between the vectors p^is and Pvis for each of the two r leptons, and 
it can be expressed in terms of the other variables. For hadronic decays of 
r's, the mmis is set to as there is only one neutrino involved in the decay. 
This reduces the number of unknowns. 

The number of unknowns (from 6 to 8, depending on the number of lep- 
tonic T decays) exceeds the number of constraints. Therefore, the available 
information is not sufficient to find the exact solution. However, not all 
solutions of this under-constrained system are equally likely, and additional 
knowledge of r decay kinematics can be used to distinguish more likely so- 
lutions from less likely ones. An example of such additional information 
is the expected angular distance between the neutrino (s) and the visible 
decays products of the r lepton. Figure [2] shows the distribution for the dis- 
tance AR=^y ('q^is — VmisV + (0vis — 0rnis)^ betwecu the directions of visible 
and invisible (missing) decay productq^ for the three distinct r decay types: 
leptonic, 1-prong hadronic and 3-prong hadronic. We incorporate this ad- 
ditional knowledge of decay kinematics as probability density functions in 
a properly defined global event likelihood to provide additional constraints 
and obtain a better estimator of M^-t- 



^For simplicity, we use the AR parametrization, although a 3-dimensional angle be- 
tween the decay products might be a more natural choice. 
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Figure 2: Example of the probability distribution functions V{AR,pr) for a particular 
value of the original t lepton momentum (pr)- These functions are used in the calculation 
of the likelihood C for three cases: 1-prong r (left plot), 3-prong r (middle plot), and 
leptonic decays (right plot) of r Icptons. These distributions depend only on the decay 
type and initial momentum of the r lepton. 



We first describe the method for the better constrained case, where both 
r's decay hadronically, and then we explain how the machinery is adjusted 
for the case of leptonic decays. When both r's decay hadronically, the sys- 
tem of Eqs. Ocan be solved exactly for any point in, for example, the (0misi, 
0mis2) parameter space. For each point in that grid, the vectors Pmisi,2 are 
fully defined and, therefore, one can calculate the distance Ai?i^2 between 
the vector Pvisi.2 and the current assumed direction of Pmisi,2- To evaluate 
the probability of such decay topology, we use AR distributions similar 
to those shown in Fig. but we take into account the dependence of the 
distribution on the momentum of the initial r lepton. If the r lepton polar- 
ization is neglected, the AR distribution depends only on the r momentum 
and decay type, but not on the source of r's. Therefore, we use simu- 
lated Z/7*— j-rr events to obtain AR distributions for small bins (5 GeV/c) 
in the initial r momentum, p, in the range 10 GeV/c<p<100 GeV/c (the 
range can be extended to both smaller and larger values). Events are sim- 
ulated using Pythia jl] supplemented with the TAUOLA package ^ for r 
decays. To simplify the calculations further, we parametrize the AR distri- 
butions by fitting them with a linear combination of Gaussian and Landau 
functions. Examples of such fits are shown as solid lines in Fig. [2j The 
PT--dependence of the mean, width and relative normalization of the Gaus- 
sian and Landau is then parametrized as Po/{x + pix^) + P2 + V?.^ + Va^"^^ 
yielding fully parametrized distributions V{AR,p), which can be used to 
evaluate the probability of a particular r decay topology. To incorporate 
this information as an additional constraint, we define the logarithm of the 
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event probability (or likelihood) as follows: 



£ = - log(P(Ai?i,Ki) X V{AR2,Pri)), (4) 

where functions V are chosen according to one of the decay types. To 
determine the best estimate for the rr invariant mass in a given event, we 
produce an M^-r distribution for all scanned points in the (0misi5 0mis2) grid 
weighed by a corresponding probability, V{ARi,Pt-i) x V{AR2,Pti)- The 
most probable value of the M^-t distribution is used as the final estimator 
of Mtt for a given event. An example of a such M^-t histogram for typical 
H^TT events of each category is shown in Fig. [31 A similar procedure can 
also be used to build estimators for other kinematic variables, if desired. 




M„ (GeV/c^) 




M„ (GeV/c ) 



Figure 3: Example of the Mtt distribution filled for all grid points in one of the H^tt 
events for each of the three decay modes. An entry for each point is weighted by its £. 
Plot on the left shows Mrr for the case of the ideal detector resolution and plot on the 
right demonstrates Mrr for the same three events in the case of the realistic detector 
resolution. All distributions arc normalized to a unit area. 

For events where one or both r leptons decay leptonically, the above 
procedure is adjusted to account for the unknown value of mmis of the 
two neutrinos in each of the leptonically decaying r's in the event. In 
this case, the scan is performed in a phase space of higher dimensional- 
ity: (0misi,0mis2!''^misi), if ouly oue of the two r's decay leptonically; or 
(0misi50mis2!''^misi5''^mis2)) if both decay r's decay to leptons. As in the fully 
hadronic mode, one can unambiguously reconstruct the 4-momenta of both 
T leptons for each point on the grid and calculate the event probability 
according to Eq. |H For simplicity, we scan uniformly in the entire range 
of kinematically allowed values of mmis, but a scan performed according to 
the mmis probability distribution function obtained from simulation may 
improve the algorithm performance. 
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3.2. Performance of the MMC Technique with Ideal Detector Resolution 

To evaluate the performance of the MMC algorithm, we use inclusive 
Z/Y^TT and gg^H^TT (with Mh=115, 120, and 130 GeV/c^) events 
produced by the Pythia MC generator supplemented with the TAUOLA 
package. All events are generated for pp collisions at ^/s =1.96 TeV. How- 
ever, the algorithm performance for events produced in pp collisions at the 
LHC is expected to be very similar to that for rr events at the Tevatron. 
Unless it is otherwise noted, we select events where both visible rs have 
Pt>10 GeV/c and ^>10 GeV is calculated as a combined transverse 
momentum of all neutrinos from both r decays). The events are categorized 
according to the decay mode of each of the two r leptons (leptonic, 1-prong 
or 3-prong hadronic) , and the rr mass is reconstructed using the appropriate 
version of the algorithm. Results for if— )-rr events with Mh =115 GeV/c^ 
are shown in Fig. H] for each of the three decay categories. In all cases, 
the peak position of the reconstructed M^-t distribution is within ~2% of 
the true mass, indicating that the assumptions used in the algorithm do 
not bias the reconstructed mass. The resolution of the reconstructed Mrr-, 
defined as the RMS of the mass distribution in the (1.0 ± 0.4) x M*™*^ 
range, changes from ~8% for events with both r's decaying hadronically to 
~13% when both r's decay leptonically. The worse resolution in the lep- 
tonic modes is due to the weaker constraints on the system. The fraction of 
events where Eqs. [3] cannot be solved for any of the grid points ranges from 
~1% to 3%, which demonstrates the high reconstruction efficiency of the 
MMC algorithm. Figure H] shows comparison of the reconstructed rr mass 
in Z/'-j* and Higgs boson events with =115 and 130 GeV/c^ when both 
r's decay hadronically. 

3.3. Effects of Detector Resolution 

To evaluate the importance of detector effects on MMC performance, we 
use the same inclusive Z/^*-^tt and gg^H^rr events and smear the pi 
and momenta of the visible r decay products accordiiig to typical detector 
resolution^ at the LHC and Tevatron experiments 

0,3 HQ- We 



assume 



3% and 10% resolutions for momenta of light leptons and hadronic r-jets, 
respectively. The pi resolution for each of the two [x- and y-) components 
is taken to be ax=(Ty=a=b GeV jsj. Note that in a real experimental en- 
vironment, the mismeasurements in lepton or hadronic r-jet momenta also 



^For simplicity, we assume Gaussian detector resolutions in this study. 
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Figure 4: Left plot demonstrates the reconstructed Mtt in H^tt events with 
Mh =115 GeV/c^ for each of the three decay categories: both t's decay hadronically 
(sohd hne); one r decays leptonically and the other one hadronicaUy (dashed hue); and 
both t's decay leptonicaUy (dashed-dotted hne). Right plot shows the reconstructed mass 
in Z/7*— >-rT and H^tt events with Mh =115 and 130 GcV/c^ in the fully hadronic 
decay mode. Results are obtained in the assumption of the ideal detector resolution. 
Each distribution is normalized to a unit area. 



lead to an additional mismeasurement in fr£. This effect is properly ac- 
counted for in our studies. Angular resolutions for visible r decay products 
of typical detectors are usually accurate enough to have no noticeable effect 
on our calculations. 

We find that mismeasurements of the momentum of r lepton decay prod- 
ucts alone have little effect on the performance of the algorithm. The M^- 
peak position and resolution are nearly unaffected and the efficiency is de- 
creased by ~3-7% as a result of mismeasurements in the momenta of visible 
r decay products, which are also propagated into the The stability of 
the peak position is related to a built-in self- correcting mechanism in the al- 
gorithm, which compensates slight under(over)-estimations in the measured 
momenta of visible decay products by over(under)-estimating the missing 
momentum, thus leading to the correctly reconstructed momentum of the 
original r lepton. 

One could expect the effects of finite ^ resolution to degrade the al- 
gorithm performance. We find that, if not taken into account, a 5 GeV 
resolution in fj^ results in a 30-40% drop in reconstruction efficiency, long 
tails in the reconstructed M^,-, and a significant degradation in the Mrr 
resolution (e.g., from ~8% to ~18% in fully hadronic rr decay mode). In 
particular, a large reduction in the reconstruction efficiency occurs because 
mismeasurements in ^ break the key assumption that the neutrinos from 
the r decays are the sole source of px: in the event (see Sec. 13. II and Eqs. [3]). 
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Figure 5: Left plot illustrates the reconstructed mass in Z^tt (solid line) and 
H^TT events with Mh =115 (dashed line) and 130 GeV/c^ (dashed-dotted line) in the 
fully hadronic decay mode. Right plot demonstrates a comparison of the reconstructed 
mass in H—>-tt events with Mh =115 GeV/c^ when both r's decay hadronically (solid 
line) and when one r decays leptonically and the other one hadronically (dashed line). All 
results are obtained by using the MMC technique in events simulated with the realistic 
detector resolution. Each distribution is normalized to a unit area. 



To mitigate these effects, the implementation of the MMC technique in a 
reahstic experimental environment has to be adjusted to allow for possible 
mismeasurements in It is achieved by increasing the dimensionality 
of the parameter space in which the scanning is performed to include the 
two components of the pr resolution (for pr^ Piy)- In this case, the 
event likelihood, C, has to be augmented with the corresponding resolution 
functions: 

£ = -log{V{ARuPrl) X P(Ai?2,Pr2) X P(A_^J X ViApTy)), (5) 

where the probability functions V{ApTx) ^nd V{ApTy) are defined as: 



exp 



2a^ 



(6) 



where a is the resolution (which we take to be 5 GeV) and Apr^ ^ are 

variations of x- or y- components of In a real experimental setup, 

the pT uncertainty can be larger in a particular direction, for example, if 

there is an energetic jet. In such cases, the uncertainty in the jet energy 

measurement will increase the uncertainty in p^ in the direction of the 

jet. These effects can be accounted for by suitably defining the x- and y- 

directions on an event-by-event basis and by choosing the appropriate ax 

and ay, which will not be equal to each other in general. 
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We evaluate the performance of the modified algorithm (with the ^ 
resolution scan) using Z^tt and H^tt events smeared with the realistic 
detector resolutions as described above. Figure \5\ shows the distribution of 
the reconstructed M^- in the fully hadronic decay mode for three samples: 
inclusive Z/7*— j-rr and gg^H^rr with Mh=115 and 130 GeV/c^. Right 
plot in the same Fig. [5] demonstrates a comparison of the reconstructed mass 
in H^TT events with Mh=H5 GeV/c^ in the case when both r's decay 
hadronically (solid line) and in the case when one r decays leptonically 
and the other one hadronically (dashed line). We find that the modified 
MMC algorithm recovers almost all lost efficiency (to the level of 97-99%) 
and significantly improves the relative Mrr resolution (to the level of ~14%). 
The reconstructed mass peak position for each of the resonances is consistent 
with the corresponding true mass. We also observe that the mass resolution 
somewhat improves (at the level of 1-2%) for events with higher ^ and/or 
higher pt of visible decay products. 

3.4- Comparisons with Existing Methods 

Figure E] shows the reconstructed M^^ distributions in H^tt events 
with M/i=115 GeV/c^ obtained by using the MMC algorithm (black his- 
togram) and the collinear approximation (red line). The events are simu- 
lated with the realistic detector resolution effects as described in the pre- 
vious section. Two categories of rr events are considered: both r leptons 
decay hadronically (left plot in Fig. [6]), and one leptonic and one hadronic 
r decay (right plot in Fig. E]). The difference in normalizations of the MMC 
and collinear approximation results reflects a higher efficiency (by a factor 
of ~1.7) of the MMC method. This is because the substantial fraction of 
events have a moderate fjT or approximately back-to-back topology and are 
non-reconstructible by the collinear approximation technique. This happens 
when small mismeasurements in ^ lead to configurations for which Eqs. [2] 
have no solution. In contrast, the MMC method resolves this problem and 
has an average efficiency of 97-99%. In addition to a better resolution in 
the core of the M^r distribution, an important feature of the MMC tech- 
nique is the absence of the long tail toward higher masses present in the 
distribution obtained using the collinear approximation. This tail is asso- 
ciated with the events of approximately back-to-back topology, where the 
collinear approximation diverges as cosA0— t-I. (A0 is the angle between 
two visible r decay products in the plane transverse to the beam line.) The 
reason for this divergency is discussed in Sec. 12. 2[ The effect is illustrated 

in Fig. [71 which shows a comparison of the ratio of the reconstructed and 
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true mass as a function of cosA0 for the two methods. In contrast to the 
colhnear approximation, the absence of long tails toward large masses in the 
MMC technique presents a significant improvement for low-mass Higgs bo- 
son searches in the H^tt channel by significantly reducing a large Z^tt 
background, which would otherwise completely overwhelm the Higgs search 
region. 




Figure 6: Reconstructed mass of the rr system for gg^H^rr events with Mh = 
115 GeV/c^ simulated with realistic detector resolution effects. Results of the MMC 
technique (solid line) are compared to those based on the collinear approximation (dashed 
line). Two categories of rr events are considered: when both r leptons decay hadronically 
(left plot), and when one of the r leptons decays to e or ^ and the other r decays 
hadronically (right plot). The difference in normalizations of the MMC and collinear 
approximation results reflects a higher efficiency of the MMC method. A long tail in the 
M-TT distribution for the collinear approximation is due to the events where the two r 
leptons have approximately back-to-baek topology. 

It is also important to point out that the algorithm efficiency and the 
shapes of likelihood C distributions are expected to be different for events 
with true r leptons and those where jets are misidentified as hadronically 
decaying r leptons. This may offer an additional handle on the backgrounds 
with the misidentified r leptons, most notably M^-l-jets and QCD multi-jet 
events, and it needs to be further investigated. 



4. Performance With Data and Monte Carlo After Full Detector 
Simulation 

To illustrate the power of the proposed method using real data, we select 
a sample of clean Z/7*— )-rr events collected by the CDF experiment in 
pp collisions at a center-of-mass energy y^=1.96 TeV at the Tevatron. We 
obtain a high purity sample of rr events in the channel where one of the r 
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Missing iUlass Calcuiator. H->TT(iVlH=115 GeV/c'^ Coiiinear Approximation, H->Tr(H„=115GeV/c^) 




Figure 7: Distribution of the ratio of the reconstructed invariant mass Mrr versus cos At/), 
where A0 is the azimuthal angle between visible decay products of the two t leptons in 
H^TT events with M/j=115 GeV/c^. Results of the MMC method (left plot) are com- 
pared to those of the coiiinear approximation (right plot). Note that the new method 
performs significantly better for nearly back-to-back topology (cos A(/i— >1), which con- 
stitutes the bulk of all rr events. 



leptons decays into a light lepton (e or /i) while the other decays into one 
of the hadronic modes. The requirement of a well isolated muon or electron 
significantly reduces QCD multi-jet backgrounds in this channel. We then 
compare the observed rr invariant mass spectrum of Z/7*— )-rr events recon- 
structed using the MMC technique with results obtained using the coiiinear 
approximation. Data are also compared with predictions obtained from 
Monte Carlo (MC) simulation. Signal events and backgrounds coming from 
Z/7*— )-ee//i/i and W + jets processes are generated by PYTHIA jsj Tune 



A with CTEQ5L parton distribution functions [10|. The detector response 



is simulated with the GEANT-3 package QCD multi-jet background 
is estimated from data by using events with lepton candidates of the same 
charge. 

4.1. Data Selections 

Full description of the CDF detector is available elsewhere jol . Here we 
only briefly discuss the algorithms and selection specific to this study. Our 
choice of selection requirements is based on the two following considerations. 
First, backgrounds in the data sample of Z/7*— >rr candidate events should 
be very low to allow unambiguous demonstration of the power of the MMC 
technique. Second, event selection should not be sensitive to potential Higgs 
boson signal to avoid biases in the currently ongoing if — j-rr analysis. To 
achieve these goals, selection requirements are chosen to be extremely tight, 
thus effective for only a small fraction of Z'/7*— )-rr and if— )-rr events. In 
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fact, the signal acceptance is reduced by a factor of ~6 compared to that 
in the ongoing search for H^tt. 



We use data collected with the inclusive electron and muon triggers [12 



These select high-pr electron and muon candidates with |77|<1. In the 
offline event selection, we require at least one reconstructed electron or 
muon candidate with transverse energy (for e's) or momentum (for /x's) 
satisfying Et>20 GeV or pt>20 GeV/c, respectively; this ensures we are 
on the trigger efflciency plateau. Data used in this study correspond to an 
integrated luminosity of 5.6 fb~^. 

Hadronically decaying r leptons are identified as narrow calorimeter clus- 
ters associated with one or three charged tracks. We refer to publication 13 
for details of r identification; the selection requirements in this study differ 
as follows. We apply track isolation (I^^) defined as a scalar sum of the 
transverse momenta of all tracks in a certain range of AR = v/A(/)2 + Ar/2 
around the highest pt (seed) track of the r candidate. The seed track is 
required to have pt>10 GeV/c and /°;^''^^^^°'^^<1 GeV/c. To suppress 
events with electrons or muons misidentified as hadronic r decays, we re- 
quire E^^/(E^*^+E™)<0.9 or (^^^'^+^^^^)/ Etr-fcsP>0-5 respectively, 
where E^'^^ and E^^^ are electromagnetic and hadronic energy deposits in 
calorimeter clusters associated with r candidates. 

For electron identification we apply the same selection criteria as in 
Ref. 13| with the following exceptions. We do not apply calorimeter iso- 
lation and require the electron to have I^^^^'^'^Kl GeV/c if the hadronic 
r in the event has one track or I^^'^'^''^=0 GeV/c if the hadronic r has 
three tracks. In addition, we require hits in the silicon vertex detector 
associated with the electron track in the case of 3-prong r decays. For 
events with 1-prong r decays, the energy of the calorimeter cluster asso- 
ciated with the electron should be consistent with track momentum such 
that E^'^ /p^'^^Kl.l. Muon identification follows the procedure described in 
Ref. [l2|. We apply the track isolation, I^^^^-^'^<1 GeV/c, to all muons 
regardless the hadronic r decay mode. 

Additional selection criteria are applied on event kinematics to suppress 
backgrounds with fake electron, muon or r candidates. The lepton (electron 
or muon) and hadronic r are required to have opposite charges to reduce the 
QCD multi-jet backgrounds. We remove events where the electron candi- 
date is identified as a conversion. We also reject events consistent with those 
coming from cosmic muons. To suppress W + jets background, we require 
the lepton (e or /x) and hadronic r to be back-to-back in the x — y plane: 
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A(l)(n/e — trk, t — trk)>2.9 radians. In addition we reject events which have 
at least one jet with Ej- >10 GeV and |?7|<3.6. Jets are reconstructed in 
the calorimeter using the JETCLU cone algorithm [3] with a cone radius 
of 0.4 in the (?7,(/)) space. To reduce Z/7*— )-ee//i/i background, we only se- 
lect events with electron transverse energy 20 GeV<-E|^<40 GeV or muon 
transverse momentum 20 GeV/c<p^<40 GeV/c. Limiting the maximum 
value of the lepton transverse energy or momentum also helps to reduce 
acceptance of potential H^tt signal. We further suppress Z/7*— )-ee//i/i 
background by looking at the invariant mass of every pair of tracks and ev- 
ery pair of clusters in the electromagnetic calorimeter. We reject an event if 
at least one pair has invariant mass consistent with the mass of a Z boson: 
i.e., it is the range 66 GeV/c^<M<lll GeV/c^ in the case of tracks and 
76 GeV/c^<M<106 GeV/c^ in the case of calorimeter clusters. 

Figure E] shows kinematic distributions for the selected Z/7*— j-rr candi- 
date events demonstrating the high purity of the sample. The MC predic- 
tions for signal and background events are normalized based on the product 
of integrated luminosity and corresponding production cross sections. 

4.2. Mass Reconstruction using the MMC Technique 

The TT invariant mass in data and MC events with full detector sim- 
ulation is reconstructed with the MMC algorithm described in Sec. 13.31 
Although Z/7*— T-rr events in our data sample have no jets, a pair of r 
leptons may be accompanied by one or more jets when different event se- 
lection requirements are applied. Therefore, we describe the pr resolution 
parametrization for events with Njet=0 and Njet>0. For this purpose, we 
only count jets with Et>15 GeV and |?7|<3.6. 

For iVjet=0 events, we perform scans for the x- and y-components of 
fjj'. The corresponding resolutions of each ^ component are parametrized 
by Gaussian distributions (Eq. E]) with width auE, which is a function of 
unclustered energjH in the event: auE=Po + PiVYl ^t- We use the same 
values of po and pi as reported in the CDF publication . 

In events with Njet>0, we consider the ^ resolution in the directions 
parallel {a\\) and perpendicular {a±) to the direction of a leading jet in 

the event. We take a±=auE and (y\\ = \J crfj^ + a'^^^, where ajet is the jet 
energy resolution which is a function of the jet Et and 77. For ajet, we 



■^The unclustered energy is defined as the scalar sum of Et for all calorimeter towers 
which are not included in electron, jet or hadronic r reconstruction. 
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Figure 8: Kinematic distributions for Z / ^* ^TT^lThVv D {I = e or /i) events selected 
from the CDF data: (a) transverse momentum of the hght lepton; (b) transverse mo- 
mentum of visible decay products of the hadronically decaying t lepton, r^; (c) missing 
transverse energy, pr, in the event [l^. 

use a simplified version (assuming, Gaussian jet energy resolution) of the 
parametrization reported in Ref. |8|. If there is more than one jet in the 
event, we project aj^t for each additional jet onto axes parallel and per- 
pendicular to the leading jet direction. These projections are then added 
in quadrature to a\\ and cr_L, respectively. Finally, we perform scans for ^ 
components parallel and perpendicular to the leading jet direction. 

4-3. Reconstructed Mass Spectrum in Data 

Figure [9] shows the rr invariant mass distribution obtained with the 
MMC and coUinear approximation methods for our data sample of Z/7*— >rr 
events. The left plot shows rr mass calculated with the MMC technique 
and compares data with the sum of background and signal predictions. The 
first bin of the distribution contains events where no solution for M^-r was 
found. We note that events unreconstructed by the MMC method are pre- 
dominantly from background processes. 
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Excellent performance of the MMC technique and its advantage over the 
collinear approximation in terms of resolution and reconstruction efficiency 
is clearly demonstrated by differences in shape and normalization of the 
Mt-t distributions in the right plot of Fig. |9l To facilitate a comparison, 
the background predictions are subtracted from the M^- distributions in 
data. Events with the reconstructed mass Mt-t->160 GeV/c^ are outside 
the histogram range and are shown in the overflow bin. The fraction of 
such events is negligible (~0.3%) for the MMC method, while it is ~18% 
fot the collinear approximation. Shapes of the distributions agree well be- 
tween data and simulation, therefore we use simulated Z/7*— s-rr events to 
estimate the resolution and efficiency achieved by the MMC technique. We 
follow the definition introduced in Sec. [3] and define resolution as the RMS 
of the M^T-/M*™^ and distribution in the 0.6-1.4 range, where M^-j. is the re- 
constructed mass and M*™^ is the generated mass. We find the resolution to 
be ~16% and the reconstruction efficiency to be ~99%, in good agreement 
with the results obtained using the simplified detector simulation model in 
Sec. 13.31 In contrast, the reconstruction efficiency of the collinear approx- 
imation method is found to be ~42%. As explained in Sec. 13. 3^ events 
where the two r's are back-to-back in the x — y plane are particularly chal- 
lenging for the collinear approximation; however, such events represent a 
major fraction of the H^tt signal. A reliable rr mass reconstruction with 
the collinear approximation is possible only for a small fraction of boosted 
rr events (with smaller angles between the r's, A0(rr), and higher values 
of pr)- The MMC method does not have such limitations, thus giving a 
substantial increase in the signal acceptance. 

5. Conclusions 

The Missing Mass Calculation method is a novel experimental technique 
proposed for reconstructing the invariant mass of resonances decaying to a 
pair of r leptons. The new method provides a substantially more accurate 
reconstruction of the mass of the rr system compared to commonly used 
techniques. Its applicability to nearly all possible final state topologies 
without loss in resolution significantly improves experimental acceptance 
for future searches for resonances decaying to a pair of r leptons. The 
new method eliminates the long tail towards higher masses present in the 
frequently used collinear approximation technique, thus promising a signifi- 
cant improvement in the sensitivity of H-^tt searches at the LHC and the 
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Figure 9: Reconstructed mass of the tt system in Z/^*^TT^lThVv v (/ = e or ^) 
candidate events using the MMC and cohinear approximation techniques: (a) tt mass 
reconstructed with MMC technique, data (points) compared to the sum of background 
and signal predictions; (b) comparison of the MMC (fihed circles are data and red line is 
the signal prediction) and collinear approximation (open circles are data and blue line is 
the signal prediction) results after subtracting the corresponding background predictions. 
Unreconstructed events are shown in the first histogram bin (Mt-t-^O). Events with 
Mtt>1Q0 GcV/c^ are outside the histogram range and are shown in the overflow bin. 



Tevatron, where the main challenge is the promotion of Z— j-rr background 
events into the Higgs boson mass region. 
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